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THE  EFFECTS  OF  THERMAL  POLLUnON  ON  RIVER  ICE  CONDITIONS 

D.  A  SIMPLIFIED  METHOD  OF  CALCULATION 


by 

S.  Lawrence  Dingman  and  Andrew  Assur 


INTRODUCTION 

In  a  previous  report,  Dingman  et  al.  (1967)  developed  a  method  for  calculating  the  tempera¬ 
ture  profile  of  a  cooling  river  below  a  source  of  thermal  pollution  and  the  length  of  ice-free  reach 
which  could  be  maintained  by  such  a  source.  Computer  programs  were  used  to  calculate  heat-loss 
rates  based  on  mean  daily  values  of  meteorological  parameters  and  to  numerically  integrate  a  com¬ 
plicated  heat-loss  expression.  The  present  paper  describes  a  simplified  approach  to  the  same 
problem,  in  which  heat  loss  is  calculated  as  a  linear  function  of  the  difference  between  water  tem¬ 
perature  and  air  temperature,  so  that  the  integration  can  be  performed  analytically.  A  simplified 
but  fairly  general  procedure  for  calculating  water-air  heat-loss  rates  on  the  basis  of  air  tempera¬ 
ture,  windspeed,  solar  radiation,  and  general  atmospheric  conditions  is  also  presented. 


BASIC  EQUATION 

The  analysis  of  Dingman  et  al.  (1967)  assumed  a  rectangular  uniform  channel  carrying  steady 
flow  and  subject  to  meteorological  conditions  which  are  constant  in  space  and  time.  Heat  load 
was  also  considered  constant,  and  the  heated  water  was  assumed  to  mix  instantaneously  over  the 
entire  cross  section  (line  source).  A  heat  balance  on  a  Eulerian  volume  element  of  the  river  be¬ 
low  the  heat  source  gave  rise  to  the  expression 


X 


-P*  ACp 


0* 


where 


Ti 


Vw  ftwCp 


(1) 


(2) 


and  X  is  distance  below  the  heat  source  (cm),  p  is  mass  density  of  water  (g  cm'*),  vw  is  river 
flow  velocity  (cm  sec"'),  h  is  river  depth  (cm),  <^p  is  specific  heat  of  water  (cal  g*'  “C*'),  Tx  is 
water  temperature  at  x  (°C),  Q*  is  heat-loss  rate  (cal  sec'*  cm’* ),  is  the  natural  river  tem¬ 
perature  upstream  ftom  the  heat  source  (°C),  is  the  heat  of  pollution  (cal  sec'*),  w  is  the  river 
width  (cm),  and  is  the  temperature  of  the  river  water  at  x  >  0  (°C). 
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Q*  is  the  net  heat-loss  rate  due  to  water-atmosphere  exchanges  by  solar  and  long-wave  radi¬ 
ation,  convection,  and  evaporation,  and  is  a  complicated  Tunction  of  water  temperature,  air  tem¬ 
perature.  windspeed,  humidity,  cloud  cover,  and  other  meteorological  factors.  Here  an  approach 
common  in  heat-transfer  problems  is  used,  and  Q*  is  approximated  by  the  expression 

Q*  =  <3o  ^  -  T.)  (3) 

where  Qq  is  heat  loss  from  the  surface  when  T*  -  Ta  is  zero  (cal  sec"'  cm"’).  T^,  is  water  tem¬ 
perature,  Ta  is  air  temperature  at  the  2-m  height  (°C).  and  (j  is  a  conventional  heat-transfer  coef¬ 
ficient  (cal  sec"'  cm"’  °C"').  Substituting  eq  3  into  eq  1.  and  introducing  an  "influence  length," 
X*  (cm),  where 


*  Pvf  ''w  ^ 

X  =  - 

Q 


gives 


X 


-X  In 


Qo  <KTx  -  Ta) 
Qo  <KT^  -  TJ 


(4) 


(5) 


or 


exp(- 


Qo  *  qrTa  -  T>) 
<1?0*  9f^i  ■"  ’ 


(6) 


Equation  6  is  a  convenient  non-dimensional  expression  for  the  downstream  decrease  of  water  tem¬ 
perature  below  a  thermal  pollution  source.  It  shows  a  simple  exponential  decay  with  distance. 

To  find  the  total  length  of  the  ice-free  reach  X  (cm)  which  would  exist  below  a  thermal  pol¬ 
lution  source  the  temperature  difference  Tg  -  Tg.  where  Tg  is  the  freezing  temperature  of  water,  is 
substituted  for  T*  -  Ta  in  eq  5: 


,  ,  Qo  +  9<Ts  -  Ta 

X  In  -  . 

Qo  +  9fTj  -  Ta) 


EVALUATING  THE  HEAT-TRANSFER  COEFPiaENT 

Equation  3  states  that  tin.  heat-loss  rate  Q*  is  a  linear  function  of  Tw  -  Tg.  The  value  of 
the  approach  outlined  here  obviously  depends  on  the  extent  to  which  this  is  true.  Many  previous 
studies  of  heat  exchange  between  natural  water  bodies  and  the  atmosphere  have  used  expressions 
similar  to  eq  3  with  some  success.  For  example,  the  Canadian  Joint  Board  of  Engineers  (1926), 
cited  in  Pruden  et  ai.  (1954),  recommended 

Q*  =  46.4  (Tw  -  Ta)  (8) 

where  Q*  is  in  cal  cm*’  day*'  and  Tw  and  Tg  are  in  °C  for  a  reach  of  the  St.  Lawrence  River  and 
Pruden  et  a/.  (1954)  developed  a  somewhat  modified  form 
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Figure  1.  Total  heat-loss  rate,  neglecting  short-wave  radia¬ 
tion,  as  a  function  of  T»r  -  T*,  with  windspeeds  at  0,  5  and  10  m  sec'*. 

Q*  =  ia3.6  +  15.5 Ta  +  43.1(Tw  -  T,)  (9) 

also  for  a  reach  of  the  St.  Lawrence  River. 

These  equations,  like  many  others  of  their  type,  are  based  on  monthly  averages  of  meteor¬ 
ological  parameters  in  the  areas  for  which  they  were  developed.  The  heat-transfer  coefficients 
and  constants  in  these  equations  have  "absorbed”  the  monthly  (in  these  cases)  averages  and  nor¬ 
mal  patterns  of  variation  of  all  meteorologic  factors  in  addition  to  air  temperature  which  affect 
heat-loss  rates.  Thus  such  equations  cannot  be  applied  with  confidence  outside  the  region  for 
which  they  were  developed,  nor  can  they  be  used  to  calculate  heat-loss  rates  for  time  periods 
shorter  than  the  averaging  periods  used  to  develop  them. 

To  provide  a  more  general  basis  for  calculating  the  constants  in  eq  3,  the  computer  program 
presented  by  Dingman  et  al.  (1967,  App.  B)  was  used  to  construct  graphs  of  total  heat-loss  rate 
neglecting  solar  radiation)  Q'  as  a  function  of  T^  -  T,  over  a  range  of  windspeeds  and  under  two 
contrasting  sets  of  general  weather  conditions  (Fig.  1).  (Appendix  A  lists  the  equations  used  in 
the  cotnputer  program  in  calculating  the  individual  heat  balance  components.)  The  windspeeds 
were  0,  5,  and  10  m  sec'*  for  1)  clear  sky  with  relative  humidity  of  50%:  and  2)  complete  cloud 
cover,  cloud  height  less  than  500  m,  ami  relative  humidky  of  100%.  Water  temperature  was  taken 
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as  OC.  The  addition  of  heat  to  the  water  body  as  solar  radiation  is  not  considered  in  evaluat¬ 
ing  the  heat-transfer  coefficient  but,  as  shown  later,  must  be  taken  into  account  in  calculating  Q  . 

In  all  cases  it  is  evident  that  the  corves  of  0  vs  TJk  "  Ta  are  close  to  linear  over  the  range 
0  <  (T„  -  Ta)  $  50  and  that  the  effect  of  complete  low  cloud  cover  is,  as  expected,  to  lower  heat- 
loss  rates.  Regression  techniques  were  used  to  approximate  the  curves  of  Figure  1,  with  the  fol¬ 
lowing  results: 


Q' 

-  105.18  +  35.08(Tw- T»); 

=  0  m  sec"* 

(10) 

Clear 

Q' 

=  220.90  +  56.27(Tw  -  Tg) ; 

‘'a 

=  5  m  sec"* 

(11) 

1  Q' 

=  .336.62  +  77.47(Tw  -  Tg) ; 

Va 

=»  10  m  sec"* 

(12) 

Q' 

=  -72.85  +  37. 10(Tg,  -  Tg): 

Va 

=  0  m  sec"* 

(13) 

Cloudy 

Q' 

=  -27.44  +  59.99(Tw  -  Tg) : 

Va 

=  5  m  sec"* 

(14) 

Q‘ 

=  17.97  +  82.88(2t  -  Tg)  ^ 

Va 

•  10  Ill  sec* 

(15) 

where  Q'  is  in  cal  cm"' day** . 

In  eq  10-15  the  values  of  the  slopes  give  the  heat-transfer  coefficient  g  which  is  effective 
at  a  given  windspeed.  Since  incoming  solar  radiation  lowers  the  heat-loss  rate  by  the  amount  Qr, 
this  value  must  be  subtracted  from  the  intercept  to  give  the  value  of  Qq  for  a  particular  case.  The 
values  of  the  slopes  and  intercepts  in  eq  10-15  increase  with  increasing  Va  for  each  set  of  general 
meteorological  conditions.  Figures  2  and  3  show  that  these  relationships  are  very  close  to  linear 
and  can  be  approximated  by 


1 

[an  =  105.18  + 23.14  Vg 

(16) 

Clear  i 

1 

[/3p  -  35.08  + 4.24  Vg 

(17) 

1  a,  =  -72.85  +  9.08  tfc 

(18) 

Cloudy 

-  37.10  +  4.58  Vg 

(19) 

where  a  is  the  intercept  and  /3  is  the  slope,  and  the  subscripts  0  and  1  refer  to  clear  and  cloudy 
conditions,  respectively. 


8AMPLB  CALCULATION 

The  following  sample  calculation  of  the  length  of  an  ice-free  reach  is  carried  out  for  com¬ 
parison  with  results  from  the  computet  pros-am  presented  by  Dingman  et  aJ.  (1967.  App.  A).  A 
power  plant  which  discharges  a  heat  of  pollution  (Qp)  of  5  10*  cal  sec**  into  a  river  of  uniform 

width  (w  -  300  m),  depth  (A  >  1.75  m),  and  velocity  (v^  >  0.75  m  sec*')  is  considered.  Assuming 
complete  mixing  at  the  pollution  site,  this  heat  input  causes  a  temperature  rise  of 

Qp 

ACpW 


_ (5)  (10*) _ 

(1  g  cm'’)(75  cm  8ec'')(175  cm)(r31 110*]  cm)(l  cal  g"*  °C**) 


1.3C 
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Figure  2*  Intercepts  in  eq  10-15  as  (unctions  o(  windspeed. 
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above  the  natural  water  temperature  of  OC.  Weather  conditions  above  the  river  ate;  cloud  cover 
0:  relative  humidity  =  50%:  windspeed  =  2  m  sec*‘:  air  temperatuie  »  -25C;  beat  gain  due  to 
solar  radiation  (Qk)  =  85  cal  cm'  day'. 

Using  eq  15 

ag  =  105.18  +  23.14(2)  =  151.46  cal  cm*' day*'  -  Q' 


and 


Qq  =  Q'  -  Qr  =  66.46  cal  cm''  day*. 

Equation  17  gives 

=  35.08  +  4.24(2)  =  43.56  cal  cm’*  day"  “C"  -  q. 

For  the  calculation  of  x*  it  is  convenient  to  change  the  units  of  q  to  cal  cm'*  sec'  'C*' : 
q  =  (43.56  cal  cm-' day*  °C")(11.157B0^] day  sec") 

-  (5.04)(10")calom-*8ec-'‘‘C‘'. 


Then 


X*  .  „  (2.604)(10')  cm. 

(5.04)  (10*‘) 

The  units  of  Qq  are  now  similarly  changed  : 

Q  -  (66.46)(1.157)(l(r*)  -  (7.69)  (KT*)  cal  cm"  sec" . 
0 

Substituting  the  values  of  x*,  Q^,  q,  T,,  Tj ,  and  Ta  into  eq  7 

(7.69)(10")  +  (5.04)  (10")  (25) 


X  =  -(2.604)  (10*)  In 

X  =  -(2.604)  (10*)  In 


(7.69)  (10-*)  +  (5.04)  (10")  (26.3) 

(7.69) (10")  +  (1.26)(10") 


(7.69)(10")  +  (1.33)(10-») 
X  =  -(2.604)(10’).ln  0.9474 
X  =  (1.409)(10*)cm  -  14.09  km. 


The  computer-calculated  value  of  X  for  this  case  is  13.18  km.  Other  calculations  by  this  simpli¬ 
fied  method  also  differed  by  less  than  10%  from  the  computer-calculated  values  (see  Table  I). 
The  simplified  method  tends  to  slightly  overestimate  ice-free  reach  lengths  when  the  initial 
polluted  water  temperature  is  above  about  1C.  This  is  as  expected,  since  the  heat-loss  rates 
calculated  by  this  method  are  based  on  a  water  temperature  of  OC,  whereas  in  the  actual 
situation  (and  in  the  computer  calculation),  heat-loss  rates  are  somewhat  higher  because 
water  temperature  is  greater  than  OC  throughout  the  open  reach. 
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Table  I.  Compaiieoii  of  ice-free  teach  lencths  calculated  by 
simplified  method  with  those  calculated  by  method  of  Dingman  et  al.  (1967)- 


Trial 

i 

mat 

2 

Trial 

3 

Trial 

4 

Qp  tea]  8ec'‘) 

5  s  10* 

Uf 

3  X  irf 

8  X  lO* 

W  (HI) 

300 

60 

200 

100 

l>  fm) 

1.76 

0.6 

5 

3 

Vw  (ni  sec**) 

176 

0.6 

1 

1.5 

^nat  CC) 

0 

0 

0 

0 

®CL  day*)! 

98.56 

98.56 

270.31 

270.31 

Ct 

0 

0 

10 

10 

H  (m)t 

500 

500 

ea 

0.403 

1.43 

0.189 

4.21 

Vg  (m  eec-‘) 

2 

6 

4 

8 

Ta  (°C) 

-25 

-10 

1 

o 

-5 

X  (km) 

(Dingman  et  at.,  1967) 

13.18 

15.44 

6.21 

186.6 

X  (km) 

(Simplified  method) 

14.09 

15.97 

6.16 

202.5 

t  See  Appendix  A  for  definitions. 


DISCUSSION 

Using  eq  10-19,  calculations  can  nowbe  made  of  heat-loss  rates,  temperature  profiles  below  a 
thermal  pollution  source  in  a  cooling  river,  and  the  length  of  ice-free  reach  maintained  in  a  river  by 
a  thermal  pollution  source,  under  fairly  general  conditions  and  without  recourse  to  a  computer. 
Equations  10-15,  and  hence  eq  16-19,  are  derived  considering  that  Tw=  OC.  This  is  true  only  at 
the  downstream  end  of  an  ice-free  reach,  but  little  error  is  introduced  by  this  assumption  unless 
Ti,  -  Tj  is  small  relative  to  T,.  Use  of  heat-loss  equations  other  than  those  listed  in  Appendix  A 
would  change  the  constants  in  eq  10-19;  liowever,  the  equations  in  Appendix  A  have  proved  ade¬ 
quate  in  several  field  tests,  and  Dingman  et  al.  (1967)  found  good  agreement  between  observed 
lengtiis  of  ice-free  reaches  and  those  calculated  using  the  equations  in  the  one  case  they  tested. 
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APPENDIX  A:  SYMBOLS  AND  EQUATIONS  USED  IN 

CALCULATING  HEAT-LOSS  SATES 

Symbols 

Symbol 

Meaning 

Units 

C 

Cloud  cover 

dimensionless  (0< 

H 

Cloud  height 

m 

Qa 

Heat  added  to  river  by  long-wave  radiation 

cal  cm"‘ 

day** 

Qar 

Long-wave  radiation  reflected  at  river  surface 

cal  cm"' 

day** 

9b 

Net  heat  lost  from  river  as  long-wave  radiation 

cal  cm"’ 

day* 

<3bs 

Long-wave  back  radiation  ftCHn  river  svface 

cal  cm*’ 

day* 

Incoming  clear-sky  solar  radiation 

cal  cm*’ 

day** 

Heat  lost  from  river  by  evaporation 

cal  cm*’ 

day* 

Heat  lost  from  river  by  convection  to  atmosphere 

cal  cm*’ 

day* 

«R 

Heat  added  to  river  by  solar  radiation 

cal  cm*’ 

day* 

®RI 

Solar  radiation  incident  at  river  surface 

cal  cm*’ 

day* 

^RR 

Solar  radiation  reflected  at  river  surface 

cal  cm*’ 

day* 

Ta 

Air  temperature  (2  m  height) 

°C 

Tw 

Water  surface  temperature 

“C 

a 

(See  eq  AS) 

b 

(See  eq  A9) 

ea 

V^>or  pressure  of  air  (2  m  height) 

rob 

®8W 

Saturation  vapor  pressure  at  temperature 

mb 

(See  eq  A13) 

<7 

Stefan-Boltzmann  constant  (l.niXlO*’) 

cal  cm*’ 

day*  °K** 

Eviitlow 

(See  Dingman  et  al..  1967  for  discussion.) 

Sidarradiatim 

“  ®RI  “  *^RR 

vidteie 

Cr,  =  CcJ0.35+ 0.061(10-01  (List.  1963)  (A2) 

=  0. 108Q„  -  (6.766) (lO**)  Q®,  (Dingman  et  al. .  1967)  (A3) 
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Lot  w«fe  n4Utfa« 

“  Ob.  -  9.  +  (A4) 

where 

=  0.970  <rTw*  (A5) 

=  (0.68  +  0.036  \/^)o  T**  (clear  sky)  (Anderson,  1954)  (A6) 

Qt  =  (a  +  5e»)oTa*  (cloudy  sky)  (Anderson,  1954)  (A?) 

a  =  0.740  +  0.025C  exp  [(-1.92)  (10‘*)Hl  (Anderson,  1954)  (AS) 

b  -  (4.9)(10-*)-(5.4)(10-*)C  expl(-1.97)(10**)H]  (Anderson,  1954)  (A9) 

Q„  >  0.03  Qa-  (AlO) 

Coivecttm  aad  evapotatloa 

Qh  =  (^3  3.9va)(Tw  ~  Tt)  (Rimsha  and  Donchenko,  1957)  (All) 

=  (l.SOko  +  6.08va)(eaw-Ca)  (Rimsha  and  Donchenko,  1957)  (A12) 

where 

kg  >  8.0  +  O.SoCTw  -  Ta)  (Rimsha  and  Donchenko,  1957).  (A13) 
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In  part  I  of  this  study,  an  expression  relating  distance  x  below  a  thermal  pollution  source  to  water  temperature  T« 
under  certain  conditions  was  developed: 

^2 

X  =  C  f  dT^/Q'’ 

n 

where  C  is  a  constant,  T|  and  Tx  are  the  water  temperatures  at  x  =  0  and  x,  respectively,  and  Q*  is  the  beat-loss 
rate.  Q*  was  represented  therein  by  a  complicated  expression  requiring  numerical  inte^ation  of  tbe  above  relation. 
Here,  the  mote  complete  equation  is  shown  to  be  well  approximated  by  Q*  =  Qo+  qli;,  -  Ta),  where  Ta  is  aittem- 
perature,  Qq  is  beat-loss  rate  when  T*  -  Ti  -  0.  and  q  is  a  heat-transfer  coefficient.  It  is  further  shown  that  for  a 
given  set  of  general  meteorologic  conditioos,  q  is  a  linear  function  of  windspeed  Va  and  is  a  linear  function  of  Va 
and  net  incoming  solar  radiation  Thus  by  specifying  general  meteorologic  conditions  v,,  and  T«  -  Ta ,  one 
can  evaluate  tbe  simple  linear  relation  for  Q*,  and  the  expression  relating  x  and  T«  can  be  easily  determined  by 
analytical  integration.  This  relation  takes  the  form  of  an  exponential-decay  curve.  Resulte  using  these  approxima¬ 
tions  are  within  10%  of  those  calculated  by  numerical  integration  for  the  cases  tested. 
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